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ABSTRACT 



Most collimated supersonic protostellar jets show a collimated wiggling, 
and knotty structure (e.g., the Haro 6-5B jet) and frequently reveal a long gap 
between this structure and the terminal bow shock. In a few cases, there is no 
evidence of such a terminal feature. We present three-dimensional smoothed 
particle hydrodynamical simulations which suggest that this morphology may 
be due to the interaction of the propagating cooling jet with a non-homogeneous 
ambient medium. In regions where the ambient gas has an increasing density 
(and pressure) gradient, we find that it tends to compress the cold, low-pressure 
cocoon of shocked material that surrounds the beam, destroy the bow shock-like 
structure at the head, and enhance beam focusing, wiggling, and internal 
traveling shocks. In ambient regions of decreasing density (and pressure), the 
flow widens and relaxes, becoming very faint. This could explain "invisible" 
segments in systems like the Haro 6-5B jet. The bow shock in these cases could 
be a relic of an earlier outflow episode, as previously suggested, or the place 
where the jet reappears after striking a denser portion of the ambient medium. 



Subject headings: Hydrodynamics - stars: pre-main-sequence stars: mass loss - 
ISM: jets and outflows 
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1. Introduction 

Protostellar jets propagate into a complex ambient medium. They are immersed in 
molecular clouds which are generally inhomogeneous and formed of many small, dense 
clouds which are sites of star formation (e.g., Bally and Devine 1994). Previous three- 
dimensional (3-D) numerical modeling of continuous and intermittent jets emanating from 
protostars were performed with beams propagating into homogeneous ambient media (e.g., 
Gouveia Dal Pino and Benz 1993, 1994, hereafter GB93, GB94; Chernin et al. 1994; 
hereafter CMGB; Stone and Norman 1994). These simulations have reproduced the basic 
features typically observed in protostellar jets: the chain of more or less regularly spaced 
knots along the beam, produced mainly in jets with intermittent ejection; or the formation 
of bow shock-like structures at the jet head, with a dense shell composed of filaments and 
blobs developed from Rayleigh- Taylor and global thermal instabihties. We here report 3-D 
numerical results for jets propagating through ambient gas with power-law density (and 
pressure) distributions, and examine their effects on the structure of the jets. 

In recent work, concerned with environmental effects on extragalactic jets, Wiita, 
Rosen and Norman (1990) and Hardee et al. (1992) have performed 2-D hydrodynamical 
simulations of adiabatic, light jets propagating into a stratified intergalactic medium. In 
contrast to such extragalactic flows, the coohng distance of post-shock protostellar jet gas is 
smaller than the jet radius and hence the assumption of an adiabatic gas is inappropriate. 
We here analyse the effects of a stratified ambient medium on radiatively cooling, initially 
heavy (denser than their surroundings) jets - a picture generally believed to be more 
appropriate to protostellar jets (e.g., Mundt et al. 1990). A more detailed description of the 
results of this investigation, with simulations covering a more extensive range of parameters, 
is presented elsewhere (Gouveia Dal Pino and Birkinshaw 1995, hereafter GB95). 
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2. 



Numerical Method 



We solve the hydrodynamics equations using a modified version of the three-dimensional 
Cartesian, gridless, Lagrangian smoothed particle hydrodynamics (SPH) code described in 
GB93, GB94, CMGB, and GB95. The ambient gas is represented by a 3-D rectangular 
box filled with particles. A coUimated, supersonic jet of radius Rj is continuously injected 
into the bottom of the box, which has dimensions of ~ 20Rj in the x-axis and 12 Rj in the 
transverse directions (y and z). The boundaries of the box arc assumed to be continuous. 
The SPH particles are smoothed by a spherically symmetric kernel function of width h. 
The initial values of h were chosen to be 0.5Rj and 0.25Rj for the ambient and jet particles, 
respectively. The jet and the ambient gas are treated as a single ionized fluid with an 
ideal gas equation of state and a ratio of specific heats 7 = 5/3. Radiative cooling (due 
to coUisional excitation and recombination) is implicitly calculated using the local, time- 
independent cooling function for a gas of cosmic abundances coohng from T ~ 10^ to 10^ K 
(Katz 1989). 

We assume an initial isothermal ambient medium (T^ = 10^ K) with density (and 
pressure) distribution stratified along the jet axis. 



where is the ambient number density, Xg is the value of x at the jet inlet (in units of 
the jet radius Rj), a is a parameter of the model and the exponent P = ±5/3 for positive 
and negative density gradients, respectively. Such profiles are consistent with the observed 
density distributions of clouds which involve protostars (e.g.. Puller and Myers 1992). A 
negative density gradient may represent, for example, the atmosphere that a jet encounters 
as it propagates through a portion of the cloud that envelops the parent source. A positive 
density gradient will occur when a jet enters an external cloud. In the simulations the 
atmosphere was held steady by an appropriate gravitational potential. 




(1) 



- 5 - 



A supersonic jet propagating into a stationary ambient gas will develop a shock pattern 
at its head, where the impacted ambient material is accelerated by a forward bow shock 
and the jet material is decelerated at the jet shock. Our models of hydrodynamical jets are 
characterized by dimensionless parameters 77 = nj{xo)/na{xo), the ratio between the jet 
and ambient number densities at Xg', Ma — Vj/ca the ambient Mach number, where Vj is the 
jet injection velocity and is the sound speed in the ambient material; and qi,s = dcooi/Rj, 
the ratio of the cooling length in the postshock gas behind the bow shock to the jet radius 
(which, for shock velocities > 90 km/s is related to the ratio measured in the postshock gas 
behind the jet shock, qjs, via qjs ~ QbsV^^ (cf., GB93). The parameters of the simulations 
were chosen to resemble the conditions found in protostellar jets. Based on observations, we 
adopt the values 77 = 1 — 10; Vj ~ 400 km/s; = 24, and Rj — 2 x 10^^ cm. We assume 
a — 0.5 and /3 = ±5/3 in eq. 1. (See GB95 for a wider range of parameters.) 

3. The Simulations 

Fig. 1 depicts the central density contours of three distinct jets when they have 
propagated a distance ~ 20Rj. All flows were initially in pressure equilibrium at the jet 
inlet, with r]{xo) = 3, na{xo) = 200 cm-^ Rj ^ 2 x 10^^ cm, Vj = 398 km s-^ M„ = 24, 
Qbs — 8.1, Qjs — 0.3, and Rj/ca — 38.2 yrs. In a) the ambient medium has positive density 
gradient (/5 = 5/3, a = 0.5), in b) it is homogeneous, and in c) it has negative density 
gradient {(3 — —5/3, a = 0.5). The cooling length parameters imply that the ambient 
shocked gas is initially almost adiabatic whereas the shocked jet material is highly radiative 
(e.g., Blondin, Fryxell and Konigl 1990, hereafter BFK; GB93). The jets reach the end of 
the computation domain at (in units of Rj/ca): a) t= 1.65, b) t=1.15, and c) t=0.95. In 
the homogeneous medium case (Fig. lb), we can identify the same basic features detected 
in previous work (e.g., GB93, GB94, CMGB). The working surface, at the head of the jet. 
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develops a dense shell (with density Ugh ~ 469n(j) of shocked jet material which becomes 
Rayleigh- Taylor unstable (GB93) and breaks into blobs. A low density cocoon containing 
post-shock jet gas is deposited around the jet beam. A dense shroud of shocked ambient 
gas envelops both the beam and the cocoon. 

In the increasing density medium (Fig. la), the cocoon is compressed and pushed 
backwards by the increasing ram pressure of the ambient gas. The bow shock shape seen 
at the head of the jet propagating into a uniform medium (Fig. lb), is replaced in Fig. la 
by an elongated structure, and the beam and shell irish/na^Xo) — 783) are highly focused. 
The high pressure of the surrounding medium promotes Kelvin-Helmholtz instabilities 
(e.g. Birkinshaw 1991) with the excitation of pinch and helical modes which coUimate the 
beam, drive internal shocks, and cause some jet twisting and flapping. A similar structure 
has been found for adiabatic jets propagating into a homogeneous ambient medium (e.g., 
GB93), but there the coUimation of the jet is a result of the hot, high- pressure, post-shock 
jet gas in the cocoon immediately behind the working surface (cf., BFK, GB93). The 
internal shocks in the beam have a typical separation of few jet diameters, in agreement 
with observations of internal knots in protostellar jets. 

In the decreasing density medium (Fig. Ic) the cocoon that surrounds the beam 
becomes broad and relaxed, the jet shows no internal shocks, and we clearly see an 
acceleration of the jet as it propagates due to the decreasing ambient density and ram 
pressure. Compared to the jet propagating into the homogeneous medium (Fig. lb), the 
beam is much less coUimated close to the head. An analysis of the density and velocity 
maps indicates a ratio of the jet to head radius e^/^ = Rj/Rh — 0.8. This decoUimation 
is significantly larger for lower-density jets (GB95). The cold shell at the working surface 
is thin, does not fragment, and has a smaller overdensity {nsh/na{xo) — 33) than the 
jet propagating into the homogeneous ambient medium, which indicates the faintness of 
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radiative shocks at the head. For an ambient medium with negative density gradient, the 
coohng length behind the bow shock (g^s) should increase as the jet propagates downstream 
and the cooling length behind the jet shock should decrease (e.g., GB93, GB95). Thus the 
bow-shock becomes more adiabatic as the beam propagates, while the jet shock becomes 
more radiative but fainter since the jet shock velocity, Vjg ~ Vj — v^s, decreases as the bow 
shock speed v^g increases. The opposite is true in the case of the ambient medium with 
positive density gradient. 

Fig. 2 depicts jets with different values of f] propagating into an a = 0.5 ambient 
medium with positive (/? = 5/3) density (and pressure) stratification. The jets reach the 
end of the compution domain at (in units of Rj/ca): a) t= 1.85, b) t=1.65, and c) t=1.25. 
The highest-density jet reaches the end of the domain earher and is least affected by the 
increasing ambient density and pressure gradients. 



4. Discussion and Conclusions 

What do these results have to do with the observed protostellar jets? Our simulations 
show that jets propagating into portions of an ambient medium with negative density 
(pressure) gradient develop a broad and relaxed cocoon and the beam rapidly accelerates 
due to the drop of the ambient density and ram pressure. No internal knots are formed in 
this case and the radiative shocks at the head are faint and will provide little radiation 
(Fig. Ic). This is consistent with the "invisible" portions of the observed outflows, e.g., the 
gaps close to the source of Haro 6-5B (FS Tau B) (Mundt, Ray and Raga 1991, hereafter 
MRR), AS 353A (e.g., Hartigan, Mundt and Stocke 1986), and 1548 C27 (e.g., Mundt et 
al. 1987). On the other hand, jets propagating into portions of an ambient medium with 
increasing density (pressure) have their cocoon/shroud compressed and pushed backwards 
by the ambient ram pressure and the beam is highly coUimated. The compressing medium 
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promotes Kelvin- Helmholtz (K-H) instabilities which cause beam pinching, twisting and 
flapping (Figs, la and 2). This favors the formation of regularly-spaced traveling internal 
knots, and an increase in the jet confinement, as required by observations, but the bow 
shock-like structure at the head disappears. 

Observed jets are, in general, well coUimated and knotty and many show one or more 
bright bow shock-like structures at the head (e.g., Reipurth 1989). In GB94 (see also, e.g., 
Hartigan and Raymond 1992; Raga and Kofman 1992; Stone and Norman 1993; Massagha 
et al. 1995), wc showed that the internal traveling knots can be a product of rapid variations 
in the jet ejection speed. We also have demonstrated that bow shocks can be formed in jets 
with continuous ejection (GB93) or with multiple outflow episodes of long period (GB94) 
and can keep their structure approximately stable in portions of the jet where the ambient 
medium is more homogeneous (GB93). There are, however, some observed jets whose 
visible structures resemble remarkably the simulated jets of Fig. 2, which propagate into 
ambient gas with increasing density gradient, as for example, Haro 6-5B-northern and HH 
24G jets (e.g., Mundt et al. 1990, MRR). 

Investigating a sample of 15 protostellar jets including those above, MRR found that 
the flows are, in general, poorly collimated close to the source (within distances ~ few 
10~^ pc) and become highly collimated further out, on scales of the order of the jet lengths 
(~ 10~^ — 10~^ pc), showing a decrease in the jet opening angle with increasing distance. 
This clearly suggests that a large scale coUimation mechanism, probably due to the external 
environment is at work (MRR). The Haro 6-5B and HH 24G jets, in particular, narrow 
with increasing distance from the source, hke the jets of Figs. 2a and b. 

The Haro 6-5B northern jet extends over Lj ~ 0.1 pc and its source is surrounded 
by an extended reflection nebula whose brightest part lies on the jet axis at ~ 7 x 10~^ 
pc. Within 10~^ pc of the source, the jet has a strong lateral expansion which is followed 
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by a gap (between ~ 1 and 1.7 x 10~^ pc) where the jet is not seen. Further out, the 
jet is recoUimated and narrows over most of its length. Its elongated structure is knotty 
and wiggling. After that, the jet becomes again invisible and then shows up in a distant 
bow-shaped nebula (~ 0.1 pc from the source). As suggested by MRR, the initial expansion 
is likely to be a signature of a highly overpressured and poorly coUimated flow. The 
appearance of the jet outside 1.7 x 10^'^ pc, where the expansion is slow {AR/Ax ~ 0.025), 
is reminiscent of the jet brightening and coUimation that we see in simulations of cooling 
jets in an ambient medium of increasing density (Figs. 2a and b). K-H instabilities here 
would drive the formation of internal shocks and small amplitude wiggles, so that the 
observed bright knots might be oblique internal shocks excited by the K-H modes (as seen 
in Fig. 2) instead of being due to intermittent jet ejection (e.g., GB94). Outside this part 
of the flow, the disappearance of the jet might arise from passage through an ambient 
medium with decreasing density (and pressure), which would lead to a low emissivity (as 
in Fig. Ic). The distant bow-shaped nebula could mark an impact between the jet and a 
last dense part of the ambient medium (e.g. GB93), or it could be a rehc of a much earlier 
outflow episode in a more homogeneous ambient medium (e.g., GB94). The morphology of 
the Haro 6-5B counterjet is radically different from that of the jet, suggesting that intrinsic 
asymmetries due to environmental effects are at work (MRR). 

Since we did not choose parameters for our simulations to fit the particular case of 
the Haro 6-5B jet, our results can provide only qualitative descriptions of the observations. 
Nonetheless, we can estimate the changes in the flow parameters that would be required to 
produce the observed expansion and collimation of the Haro 6-5B jet. The coUimated part 
of the jet flow has and expansion rate [AR/ Ax)obs — 0.025 (MRR), and in our simulations 
of the 77 = 3 jet propagating into an increasing pressure environment (Fig. la. Fig. 2b) 
we find AR/Ax ~ 2x smaller. This suggests that a reduction of the jet overpressure by a 
factor ~ 4 would cause a better match between our simulations and the observed flow. We 
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note that the ambient temperature used in our simulations (and thus Pa) is somewhat larger 
than the observed values in star forming regions, an approximation that has been used 
in previous work (e.g., BFK, GB93, GB94) partially to mimic the presence of a magnetic 
pressure. A decrease of in our simulations would be consistent, for example, with the 
required decrease of Pa in order to fit the observed expansion and narrowing of the Haro 
6-5B jet. 

Large amphtude side-to-side wiggles and knots are also observed in HH30, HH83, 
HH84, HH85, HHllO, and other jet systems (e.g., MRR, Reipurth 1989). Such features 
could also be due to the mechanism discussed here. We note, however, that some of 
these jets keep broadening slowly with distance (e.g., HH30) or have an approximately 
constant width over a section of their length (e.g., HH83). This would suggest less efficient 
confinement by the ambient medium in these cases, perhaps because of a larger initial 
density ratio between jet and ambient gas, as in the high rj jet of Fig. 2c, or a smoother 
density stratification of the ambient medium (smaller a, see eq. 1). 

Although, on average, the jets of Fig. 2 are decelerated by the increasing ram 
pressure of the ambient medium, we find that the velocity of advance of their working 
surfaces, Vws, has an oscillating pattern. Since the density ratio f] decreases with distance, 
Vyjs ^ + {r]e)~^^'^] (see, e.g., eq. 1 of GB93) initially decreases as the jet propagates 

downstream with a ratio of jet to head square radius e = {Rj/RhY ~ 1. Later, as the jet 
head is compressed by the ambient medium, the increase of e more than compensates for the 
decrease of 77, and the jet head accelerates as it narrows. Eventually, the further decrease 
of 77 dominates the equation and the working surface decelerates again. For example, the 
r]{xo) — 3 jet of Fig. 2b follows the Vy,s{^) relation above within ~ 30%. The internal 
shocks of Fig. 2 are regularly spaced (at 1 to 2 Rj) and travel downstream with a similar 
variable velocity pattern. Their speed gets closer to v^jg as they approach the jet head. 
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Recent measurements of the tangential velocity of the knots of Haro 6-5B (Eisloffel 1993) 
indicate an oscillating pattern in qualitative agreement with the results above. Further 
measurements of variations of the proper motions of knots in systems like HH24G or HH83 
would be an interesting check on our predictions. It is also important to note that the 
variations of the shock velocity at the working surface imply variations in the intensity of 
emission from the associated cooling regions behind the shocks. 

We are indebted to the referees (Colin Masson and an anomymous referee) and to 
Tom Ray for their fruitful suggestions. Also EMGDP is thankful to R. Mundt for his kind 
hospitality during the Workshop on Disk and Outflows from Young Stars at the Max-Planck 
Institute in Heidelberg and his relevant and fruitful comments on this work. The simulations 
were performed on HP Apollo 9000/720 and DEC 3000/600 AXP Workstations whose 
purchase was made possible by the Brazilian Agency FAPESP. This work has been partially 
supported by the Brazilian Agency CNPq. 
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Figure Captions 

Figure 1. Central density contours of jets with rj{xo) — 3 and Ma — 24. In a) the 
ambient medium has positive density gradient (/3 = 5/3, a — 0.5), in b) it is homogeneous, 
and in c) it has negative density gradient (/3 = —5/3, a = 0.5). The coordinates x and z 
are in units of Rj. The contour hues are separated by a factor of 1.3 and the peak of the 
density scale is 800/na{xo). 

Figure 2. Jets with distinct r}{xo) in an a = 0.5 ambient medium with positive 
stratification, a) r]{xo) — 1, b) 77(^0) = 3, and c) 77(^0) = 10. The initial conditions are the 
same as in Fig. 1. The maximum density of the shells at the jet heads are nsh/n-aixo) ~211 
(a), 638 (b), and 2540 (c). 



